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o ■ Abstract 

i-> : 

r"! , In the last few years a number of works reported the appearance of thick branes with 

internal structure, induced by the parameter which controls the interaction between 
two scalar fields coupled to gravity in (4,1) dimensions in warped space-time with one 

^ I extra dimension. Here we show that one can implement the control over the brane 

thickness without needing to change the potential parameter. On the contrary, this is 
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lO ■ going to be done by means of the variation of a parameter associated to the domain wall 
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degeneracy. We also report the existence of novel and qualitatively different solutions 



t^^ ' for a critical value of the degeneracy parameter, which could be called critical Bloch 

QQ ■ branes. 
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1 Introduction 

The problem of the so called thick branes have received a considerable amount of attention 
along the last years [l]-[23]. As observed by Campos [I] a few years ago, some kinds of 
two interacting scalar fields potentials can be used in order to describe the splitting of thick 
branes due to a first-order phase transition in a warped geometry. In that work, Campos 
discussed the effect by studying a model without supersymmetry. In a recent work, Bazeia 
and Gomes j2] discussed the appearance of thick branes by using a two interacting scalar 
fields model that can be naturally incorporated in supergravity, what was also done by Eto 
and Sakai [5]. In our work, by using the very same model discussed by Bazeia and Gomes, 
we show that the reported Bloch branes [2] have, in fact, more general soliton solutions and, 
as a consequence, the resulting brane can have a much richer structure as, for instance, a 
degeneracy controlling parameter [3l [24|. It is shown that for a convenient choice of this 
parameter a double wall structure and a corresponding thicker Bloch Brane, what we call 
degenerate Bloch branes shows up. It is important to remark that both in the case of 
Campos as in the Bazeia and Gomes, the splitting of the thick branes is controlled by the 
potential parameters. Here, instead, the splitting is controlled by means of a parameter 
which is not present in the Lagrangian density, on the contrary it appears in the solutions 
as a shape controlling parameter. In fact, as asserted in above, there is a degeneracy in the 
solution because, in spite of the value of this parameter the energy of the field configuration 
is precisely the same. In view of this we call it a degeneracy parameter [21]. In fact, 
Bazeia and collaborators [25] have introduced a one scalar field model with similar properties 
like the existence of double- walls solutions and thick branes [6] but, once again, the brane 
thickness control is done through a fine tuning of the potential self-interaction parameters 
themselves. Here the idea is to get a more robust way of controlling the thickness, and the 
correspondent distance between the walls whereas preserving a supersymmetric structure. 
Finally we introduce a special solution, at the critical value of the degeneracy parameter, 
which presents a quite different and interesting behavior for the brane. 

On the other hand finding exact classical solutions, particularly solitons, is one of the 
problems on nonlinear models with interacting fields [2n]-|H]. As pointed out by R. Ra- 
jaraman and E. Weinberg |3T], in such nonlinear models more than one time-independent 
classical solution can exist and each one of them corresponds to a different family of quan- 
tum states, which come into play when one performs a perturbation around those classical 
solutions. 

In order to deal with the systems we are going to work here, it is common to use the so 
called trial orbits method |32] , which is a very powerful one presented for finding exact soliton 
solutions for non-linear second-order differential equations of models with two interacting 
relativistic scalar fields in 1+1 dimensions, and it is model independent. A couple of years 
ago one of us presented a method for finding additional soliton solutions for those special 
cases whose soliton solutions are the BPS ones [35] and in the last year that approach was 
extended, allowing more general models [311 EZl SI]- This last approach is the one we will 
use along this manuscript. As a consequence, we present more general soliton solutions and 



how they are intrinsically related to degenerate and critical Bloch branes. Furthermore, 
once we have found more general solutions which engenders thicker branes, we discuss the 
influence of those solutions in the warp factor and in the fluctuation of the metric around 
those classical solutions. We also compare our results with those obtained in [2j. 

This work is organized as follows: In the second section we present the model we are 
going to work with and review the approach introduced in the reference [35] to find classical 
soliton solutions. In that section we also obtain the warp factor in a general form. In the 
third section, a variety of soliton solutions that have been found up to now, are constructed 
by using the method of the second section; we also present the warp factor for each set of 
solutions. The following section is devoted to discuss the stability and the zero modes for 
each set of soliton solutions in the context of the brane worlds scenario. Finally, we address 
final comments on the soliton solutions and their consequences and applicability in the brane 
world scenario. 

2 Gravity coupled to two interacting scalar fields: an- 
alytical solutions 



The action we are going to work with, is the one of a five-dimensional gravity coupled to 
two interacting real scalar fields, which can be represented by 
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where g = det {gab) and for granting that the four-dimensional space has the Poincare in- 
variance, it is usually used that 



ds'^ = gab dx"'dx = e^ ^^'rf^ydx^dx^ — dr'^] a, 6 = 0, ..., 4, 



(2) 



where r is the extra dimension, rf^i, the usual Minkowski metric and e^^*-'"-' is the so called 
warp factor. A usual hypothesis is that the warp factor depends only on the extra dimension 
r. Besides, one can also assume that the scalar fields depend only of the extra dimension r. 
Under these assumptions, one can determine the resulting equations of motion for the above 
system as [21 [131 [Tl] 
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As can be demonstrated [I3l[2], the above set of second-order non- linear coupled equations 
has another set of first-order differential equations which shares solutions with it, and that 
is given by 

d<P _ dW{(t>.x ) dx_dW{(t>,x) dA 2, 
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provided that the potential V {(f), x) is restricted to be of a given class of potentials which 
can be written in terms of a kind of superpotential as 



V{(^,X) 



dw{4>,x) 



+ 



dW{(^,x) 
dx 



lw{<p,x? 



(6) 



Note that, once one has solutions of the first two first-order equations for the interacting 
scalar fields, it becomes a simple task of integration to get A (r) and, as a consequence, 
determine the warp factor. Now, in order to go further on the analysis it is important to 
work with a concrete example which, when available, should be one with exact analytical 
solutions [TB] . So, we will work with a superpotential equivalent to that used by Bazeia and 
collaborators [2]. The idea is to show that there are other solutions which were not analyzed 
in reference [2], and which present quite interesting features. Particularly, we will show that 
some of them allow one to control the behavior of the warp factor without performing the 
kind of restriction over the potential parameters as the one did in |2]. The superpotential 
we are going to work with along this manuscript is 



W(4,,x)=<t' 
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which becomes equal to that considered in [2] by choosing a = 1, A = — 1 and /j, = —r. 
From now on, in order to solve the first-order differential equations, we follow the method 
of reference [35] instead of applying the usual trial orbits method |32], [38]. For this we note 
that it is possible to write the relation d(f)/W^ = dr = dx/W^, where the differential element 
dr is a kind of invariant. Thus, one is lead to 



d^^W^ 
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This is in general a nonlinear differential equation relating the scalar fields of the model. If 
one is able to solve it completely for a given model, the function (x) (in fact, it will be the 
equation for a generic orbit) can be used to eliminate one of the fields, rendering the first- 
order differential equations uncoupled and equivalent to a single one. Finally, the resulting 
uncoupled first-order nonlinear equation can be solved in general, even if numerically. By 
substituting the derivatives of the superpotential ([7]) with respect to the fields in ([8]) we have 
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which can be rewritten as a hnear differential equation 

dp A 



dx /^X 



X, (10) 



by the redefinition of the fields, p = cfP' — a?. Now, the general solutions are easily obtained 
as 

p(x) = 02 - a^ = CO x'/^ - Y^^ ^'' f°^ ^ ^ 2/i, (11) 

and 

p{x) =ct)^-a^ = x^[Hx) + ci], for A = 2/i, (12) 

where cq and ci are arbitrary integration constants. We substitute the above solutions, 
for instance, in the differential equation for the x field, obtaining the following first-order 
differential equations for the field x(^) 

'^ = ±2i,xJa^ + c^X^'^-Y^^\ A ^2/., (13) 

and 

^ = ±2/xxVa2 + x2[ln(x) + ci], A = 2/i. (14) 

As a matter of fact, in general, an explicit solution for each one of the above equations 
can not be obtained, but one can verify numerically that the solutions belong to the same 
classes, and some of those classes of solutions can be written in terms of analytical elementary 
functions. In those last cases one is able to obtain the several types of soliton solutions we 
discuss in the next section. 

Now, let us discuss a bit on the form of the warp factor e'^^^^\ for which it is necessary to 
compute the function A (r) by integrating the second equation in (5). Once we will consider 
many situations, some of them with expressions much more involved than those studied in 
[2] , it should be very convenient to express it in a general form in terms of the field itself, 
instead of a function of the spatial variable. Furthermore, that would allows one to make 
qualitative considerations about the behavior of the warp factor. In order to put this idea 
in a concrete form, we will use the orbit equations (TTTl [T2l) and manipulate the equation for 
A (r) in order to write it in terms of the field x (^)- Fo^' this we start by noting that after 
eliminating the dependence of A (0, x) by using the orbit equation (0 = (x)), one obtains 

dA (r) _ dA (x) dx _ dA (x) dW (0 (x) , x) _ 2 
dr dx dr dx % ^'^ i<PiX) ,X) , U^J 

which leads to 
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where W^ (0 (x) > x) — — g • Now, substituting the superpotential of the model under 
analysis we get 



dAix) 1 



^ I 0(X) 2 1, 2 



dx 3 fix 

which after some simple manipulations using the orbit equation conduces to 
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Finally we can perform the integration over the field x? obtaining 

A (,, = „. + (!i^) ,„(,, _ I (A^) ,^ _ (I) ,«). for A ^ 2,. (20) 

and 

^ (X) = «i+ (^) In (X)-^^^^^I^ X'--^ X' (in (x) " ^ , for A = 2^, (21) 

where cto and ai are arbitrary integration constants, which are going to be chosen to ensure 
that A{r = 0) = 0. It is important to remark that the above solutions are completely general 
for this model and, as a consequence, can be used to get the warp factor for an arbitrary 
choice of the potential parameters, even for the cases where its solution can not be obtained 
through analytical elementary functions. The above general approach can be checked, for 
instance with case studied by Bazeia and Gomes [2], for which 



x(l/) = ±\/^-2sech(2sy), (22) 



and, from above, one obtains 



A{y) = — Ul-3 s)] tanh (2 syf -2 In (cosh (2 s y)) , (23) 
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where we have used the original variables and parameters defined in PJ. This is precisely 
the result obtained in that work through direct integration in the spatial variable. 

From the expressions obtained above, one can clearly see that the behavior of the warp 
factor is very sensitive to the one of the field x- For instance, when x {'^) changes very slowly 
in a given region, so it will happens with A (r). In certain manner, one can guess the behavior 
of the warp factor, simply by observing that for x- 



3 Soliton solutions and their warp factors 

In this section we explore in some extent the solutions for the equation (fT3|) . This is done 
by presenting those resulting soliton solutions for the model under consideration and also 
by obtaining explicitly the warp factor for each set of soliton solutions in the brane scenario 
under analysis. Finally we compare our results with those offered in [2]. 

Before proceeding with this program, we would like to stress that the model we are 
working with admits a particular set of solutions which can not be obtained from the method 
described in the previous section. That set of classical solutions, which could be called 
isolated solutions because is characterized by Xi{^) = 0? such that there is no sense in 
writing the differential equation ([9]) for this case and, consequently, do not furnishes any 
internal structure for the brane [2j • Even though, the system admits a soliton solution given 
by 0/(r) = ±a tanh(Aar), where the (lower) upper sign refers to a (anti-)kink solution. 
We will not consider this case in detail, once it is effectively a one field model, and we are 
primarily interested in the two fields nontrivial solutions. 

3.1 Bloch walls 

The usual set of solutions, baptized as Bloch wall in [2j, can be obtained by means of 
the method described in the previous section. It is obtained when we take Cq = in the 
expression flT3l) . In this case that equation can be solved analytically for any value of A and 
/i, provided that \ > 2fi in order to keep the solution real. In this case we get the following 
solution for x(r) 

XBwi'r) = ax sech(2 /iar). (24) 



One can observe that this solution vanishes when x -^ ±oo. The corresponding kink-like 
solution for the field 0, is given by 

4>Bw{^) = ict tanh(2 /i a r), (25) 

which connect the vacua of the potential. In this case the warp factor of the configuration 
is the one presented in the previous section in equation (l23l) . We call this type of domain 
wall as BW domain wall to distinguish it from other types of domain wall solutions we are 
going to present for this model. 

At this point it is interesting to note that, some of the solutions we are going to explore 
in the next sections are, in fact, in a different range of the potential parameters as compared 
with the ones considered in [2]. For instance, when we consider the A = /i, in terms of 
the parameters used by Bazeia and Gomes, this case would corresponds to consider r = 1. 
However the range of validity of the solutions used by them is < r < ^ 



3.2 Degenerate Bloch walls 



Others soliton solutions can be found when one considers the integration constant cq 7^ 0. 
It was found in the reference [35] that at least in three particular cases the equation flT3|l 
can be solved analytically. For cq < —2 and A = /x it was found that the solutions for the 
x{r) field are lump-like solutions, which vanishes when r — >■ ±cx3. On its turn, the field 0(r) 
exhibits a kink-like profile. 

These classical solutions can be written as 
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^ycl — 4 sinh(2 ^ar) 



a/cq — 4 cosh(2 fiar) — Co 
and its warp factor is given by 
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for A = /i, Co < -2, (26) 



for A = /x, Co < -2. (27) 



2 a ( Cq — Co a/cq — 4 cosh(2 ^lar) — 4 < 
9 ( a/cq — 4 cosh(2 fiar) — Co 



(28) 

where, as we anticipated above, A^q, will be chosen in order to get e^"^*^"-* = 1, for plotting 
convenience. 

An interesting aspect of these solutions is that, for some values of Co < —2, 0£,^vi/('") 
exhibits a double kink profile. We can speak of a formation of a double wall structure, 
extended along the space dimension. In the Fig.l we compare the case studied in [2] to 
some typical profiles of the warp factors in the case where A = /i, both when co is close to 
its critical value (co = —2 in this case) and far from it. One can observe the appearance of 
a more pronounced flat region, where one could speak of a Minkowski- type metric. In fact 
this "Minkowski sector" becomes larger as co approaches its critical value. 

Similar behavior is also noted in the classical solutions for A = 4/i and co < 1/16. In this 
case the field x{^) has a lump-like profile given by 
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with the corresponding warp factor being 
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1 + 32 a^co + -\/l — 16co cosh(4 fiar) 
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and, once more we choose N^ in order to get e^^^^^ = 1. 

In this last case 4>DBwi''^) ^^^^ presents a double kink profile for some values of cq. In 
Fig. 2 it can be seen the rising of two peaks at the extremum of the flatten region and an 
increasing of the distance between the two peaks in the warp factor as cq approaches its 
critical value (cq = 1/16 in this case). 

3.3 Critical Bloch walls 

Finally, a very interesting class of analytical soliton solutions were shown to exist when one 
takes X = fi with the critical parameter cq = —2 and for A = 4/i with the critical parameter 
Co = 1/16, in the equation flT3l) . The novelty in these cases is the fact that both, the x(^) and 
the 0(r) fields present a kink-like profile and the warp factor presents a remarkable behavior, 
what can be noted from Fig. 3, where it is evident the existence of two " Minkowski-type" 
regions, separated by a transition one. This argument is reinforced by the behavior of the 
energy densities of the soliton configurations, as can be noted in the Fig. 4, as well as from 
the stability potential (see Fig. 5 and 6). 

We call this set of solutions as CBW domain walls. For X = fj, and cq = —2 the classical 
solution for the x{^) fi^ld can be shown to be 



XcBwi^) = 2 ^^ ^ *^^^ if^^^r)] 



(32) 



and the solution for the 0(r) field is given by 
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The corresponding warp factor is 
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For Co = 1/16 and A = 4/i, the following set of domain walls is obtained 
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and 
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The warp factor, on its turn, is found to be given by 
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(37) 
X exp < — tanh(2 fxar) [l6 o^Cq tanh (2 /i a r) =l= (1 + 32 a^Co)] [■ . 

4 Stability and zero modes 

In general it is quite hard to take into account a full set of fluctuations of the metric around 
the background in a model where gravity is coupled to scalars. This happens as a consequence 
of a very intricate system of coupled non-linear second-order differential equations [13l [HI 
[T6l E]. Fortunately however, there is a sector where the metric fluctuations decouple from 
the scalars, and it comes to be the one associated to the transverse and traceless part of the 
metric fluctuation [131 [H]- This can be shown, if one introduces a metric perturbation like 



ds^ = e^^^'') [r]^^ + e h^^) dx^'dx'' - dr^, (38) 

and perform small fluctuations on the scalar fields, (p ^ (f){r) + e (j){r, x^) and X "*■ X (^) + 
e x{r,x^), with h^i, = h^^{r,x^), and e is small perturbation parameter. Now, keeping 
the terms in the action up to the second order in e, as done originally by DeWolfe and 
collaborators [13], as well as by Bazeia and Gomes in the case of two scalar fields [2], one 
gets the following set of coupled equations for the scalars fluctuations 

dr dr dr"^ dcp"^ d(j)dx 2 dr dr ' 



(39) 



dr dr dr'^ dx^ d(j)dx 'i dr dr 



and, for the metric fluctuations, one obtains 



+ In^y^'^Or (,-.„.) + le^%. (I ^ + I ,) . 0. (40) 

One can simplify this last equation by choosing a transverse and traceless h^^, which is done 
through the use of the projector 

P^lua/3 = - ij^^ia T'ul3 + T^ tiP T'ua) " -'"'^lu T^a/S, (41) 
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where tt^j, = rj^^i, ^. In other words, by using that h^^i, = P^uafs ^°^, one arrives at 

dr^ ^ dr dr " ^"^ ^ ' 

Now, performing a sequence of function redefinition 

/i^, = e^*^-e-i^^^„ (43) 

and variable transformation z = J e~^^'^^dr [2j, one can recast the above equation into a kind 
of Schrodinger equation 

- ^ + U^ff (z) ^^. = k'ij,., (44) 

where the effective potential is defined by 



9 /dAy 3d^A 
4 \lh) ^21^ 



Ueff{z) = -[ — ] +--^. (45) 



The above differential equation can be factorized as 



d ?)dA\ ( d 3(iy4\ , , ,9 



It can be shown that A;^ is positive or zero since the resulting Hamiltonian can be factorized 
as the product of two operators wich are adjoints of each other [21]. So the system is stable 
against linear classical metric fluctuations. 

Regarding the stability of the system against classical linear fluctuations of the scalar 
fields we note that the equations (15^ are quite hard to be analyzed due to the couphng of 
the fields among themselves and with the metric fluctuations (The term in the right hand 
side of the equations fl39l) ). Thus, one can try to simplify the problem by considering only 
the fluctuations of one of the scalar fields. If we consider, for instance, only the fluctua- 
tion of x('^)) we are left just with the second of the equations fl39l) . Performing the field 
transformation x(r) = e^"^'^ e~^^(^{z) we obtain 






-+ f/,^n^) + e'^^^^TrT)C = A, (47) 



where t/e// (-z) is given in the equation (H5l) and we have taken into account only the zero- 
mode of the metric fluctuation. This can be done as far as the potential f/g// (-z) , the 
one responsible to localize the gravity, supports only the zero mode as a localized state. 
Unfortunately, we have not been able to factorize the equation Wh as a product of two 
operators wich are adjoint to each other, as in equation fH6|) . If that is possible the system 
is also stable against the fluctuations of at least one of the scalar fields. Furthermore, in 
general the factor e^^'^^^ {(f'V j d}^^ is not positive for all values of the variable 2;, thus we can 
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not guarantee that the spectrum is positive semi-definite. As far as we know the question 
regarding the stabihty of the scalar fields is an open problem. In the reference (??) this 
question was thoroughly examined for the situation where only one active scalar field is 
present and all the metric fiuctuations modes were fully taken into account. The authors 
show that in some cases the stability of domain walls can be proven, that is p^ > 0, although 
the effective potential can not be factorized. Unfortunately, the case we are studying here 
does not belongs to any of those cases. 

Returning to the analysis of the metric fiuctuations equation, we remark that the zero 
mode coming from the equation (44) grants the existence of massless four-dimensional gravi- 
tons [H [ini [HI [2] . In general the shape of the zero mode is quite similar to the warp factor so 
that one can think that there is some relation between them. In the next we will show that 
these two quantities really present the same generic shape. With this in mind, we start from 
equation ( H2i) . redefine the function h^i, as hfj^^y = e"^'^'^''^^^,^ (r), and obtain the following 
equation 



dr"^ \ dr^ \dr J 

which, for the case of the zero mode (/c^ = 0) can be rewritten as 



d dA\ ( d dA\ ^ , , 

+ 2-^ -- r + 2-^KM- = 0> (49) 



dr dr ) \ dr dr ) 

and finally one can see that the zero mode solution, apart from a normalization factor, is 
precisely the warp factor 

e^ = iVoe^^('^V. (50) 

In terms of the coordinate r the effective potential which localize the gravitation in the 
brane is written as 

Obviously the above potential is equal to the one in the z variable, it will be a kind of 
re-scaled one. However the general shape and characteristics in both variables are the same. 
The stability potential is represented in figures 5-7. In Fig. 5 we compare the behavior of 
one of our degenerate cases with those of Bazeia and Gomes [2] . Fig. 6 shows clearly that 
the structure of the potential in a situation where two interactive regions are separated by 
an approximately zero force one. Finally, in Fig. 7, we see that those separated potentials 
recombine into a single one. 

The essential idea in our work is to show that the situation is much richer than that ana- 
lyzed in [2], and that from a complete set of solutions as the one we present here, important 
consequences for the warp factor structure and, consequently, for the brane world scenario 
show up. One can cite for instance the fact that one can control, by means of a parameter 
which is not present in the potential V (0, x), the region where the metric is approximately 
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flat. Furthermore, in a given critical case, there are two of these regions, separated by a 
transition one (see Fig. 3). 



5 Final remarks 

In this work we analyze the impact of a general set of solitonic solutions over the charac- 
teristics of some models presenting interaction between two scalar fields coupled to gravity 
in (4,1) dimensions in warped space-time with one extra dimension. Essentially, we explore 
a larger class of solutions of a model recently studied [2]. In doing so, we have discovered 
a number of interesting features as, for instance, a kind of type-I extreme domain wall as 
classified in [28], when we have dealt with what we have called critical domain walls (see 
Fig. 3). 

One very important consequence of our study is that the thickness of the domain walls 
can be controlled by means of an external parameter (regarding the scalar fields potential), 
and this can be done without changing the potential parameters, in contrast with is done in 
other models [21 [6] . 

Furthermore, one can observe the appearance of a controllable flat region in the warp 
factor, where one could speak of a Minkowski-type metric region (see Fig. 2 and 4). In fact, 
in Fig. 4, where the energy density is plotted, we see clearly in the case with cq = —2.00001, 
that the region where negative energy densities show up is outside of the " Minkowskian" 
one. Thus one could speculate about a possible confining mechanism for the bulk particles 
in that internal region. 
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Figure 1: Warp factor appearing in the reference [2j, with the parameters used there: r = 0.05 
(thick sohd hne), 0.1 (dashed hne) and 0.3 (thin sohd hne) (left). Warp factor for the case 
where: a = 1 A = /i, Cq 7^ — 2 and Cq = —2.00001 (thick sohd hne); —2.005 (dashed hne); 
—3.0 (thin sohd hne) (right). 
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Figure 2: Warp factor for the case where: a = 1 A = 4/i, 
(thick sohd hne); 1/17 (dashed hne); 1/200 (thin sohd hne). 
peaks, signahzing a richer structure for the zero mode. 
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Figure 3: Warp factor for the case where: A = 4/i, Cq = 1/16, /i 
hne), a = 1.2 (dashed hne) and the case where a = 2 and fj, = 0.2. 
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Figure 4: Energy density of reference [2] r = 0.05 and 0.30 (left). Energy density for the 
case where: a = 1 A = /i = 1, Cq 7^ — 2 and cq = —2.00001 (sohd hne); —4.0 (dashed hne) 
(right). 
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Figure 5: Comparison of the stabihty potential of reference [2] (left) with r = 0.05; 0.1; 0.3 
and ours with A = /i, a = 1 and Cq = —2.001 (thick sohd line); —2.1 (dashed line); —2.5 
(thin solid line). 
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Figure 6: Stability potential for the case where Cq = 1/16.0001 with a = 1 (thin dashed 
line), and a = 1.2 (thick solid line) 




Figure 7: A typical stability potential for the critical case, both when A = ;U as A = 4/x. 



